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Abstract 

Two coarse-grained models for polymer chains in dense glass-forming poly- 
mer melts are studied by computer simulation: the bond fluctuation model 
on a simple cubic lattice where a bond-length potential favors long bonds is 
treated by dynamic Monte Carlo methods, and a bead-spring model in the 
continuum with a Lennard-Jones potential between the beads is treated by 
Molecular Dynamics. While the dynamics of both models differ for short 
length scales and associated time scales, on mesoscopic spatial and temporal 
scales both models behave similarly. In particular, the mode coupling theory 
of the glass transition can be used to interpret the slowing down of the under- 
cooled polymer melt. For the off-lattice model, the approach to the critical 
point of mode coupling is both studied for constant pressure and for con- 
stant volume. The lattice model allows a test of the Gibbs-Di Marzio entropy 
theory of the glass transition, and our finding is that although the entropy 
does decrease significantly, there is no "entropy catastrophe" where the fluid 
entropy would turn negative. Finally, an outlook on confinement effects on 
the glass transition in thin film geometry is given. 

Contribution to the II Workshop on Non Equilibrium Phenomena in Supercooled 
Fluids, Glasses and Amorphous Materials, Pisa 1998, submitted to J. Phys.: 
Condens. Matter on Sept. 27, 1998 

I. INTRODUCTION 

The glass transition from the supercooled melt to the amorphous solid still is a puzzle 
. While the viscosity rj(T) increases about 15 orders of magnitude (T g is empirically 
defined by rj(T = T g ) = 10 13 Poise) in a narrow temperature interval, the change in the static 
structure factor S(q) (q is the scattering wavenumber) is rather minor, see e.g. ||. What is 
then the physical basis that many different systems show a nearly universal phenomenology 
of relaxation {e. g., increase of relaxation time r or rj(T) [t)(T) oc r] described by the Vogel- 
Fulcher-Tammann (VFT) law, lnr oc Inrj cx Eyp/(T — Tvf) 0, relaxation functions obey 
the Kohlrausch- Williams- Watts (KWW) law, <p(t) oc exp [-(t/r)^], etc. 0}? 

The answer to this question is controversial: e.g., there is the idea of an "entropy catas- 
trophe" static phase transition underlying the glass transition {the configurational entropy 



1 



of the undercooled fluid vanishes at T < T g 0, with possibly @ T = Typ}. Mode coupling 
theory 0, on the other hand, yields a (rounded ||) dynamical transition from an ergodic 
fluid (atoms easily can escape from the cages formed by their neighbors) at a "critical" tem- 
perature T c > T g to a nonergodic state ("structural arrest"). Still other concepts attribute 
the increase of r upon cooling to the increase of a glass correlation length ||, in analogy 
with spin glasses [{lOfl . 

Polymer melts can be held in very good metastable equilibrium in their supercooled 
state, and their T 9 's occur in an experimentally convenient range, and thus a wealth of 
experimental data exists. Thus it is desirable to approach these systems also by theoretical 
modeling via computer simulation JO]. However, due to the large size of these polymer 
coils (they exhibit nontrivial structure from the A scale of covalent bonds up to the gyration 
radius which can be 100 A) and the many decades of time scales spanned by the relaxation 
times of their motion simulation of polymers is very difficult, and requires to use simplified 
efficient models ||11||. Two such models will be defined and compared in the next section, 



while Sec. 3 describes the tests of mode coupling theory performed with these models. Sec. 
4 briefly describes tests of other theories and contains also some concluding remarks. 



II. COARSE-GRAINED MODELS FOR POLYMERS ON THE LATTICE AND IN 

THE CONTINUUM 

The basic idea is to eliminate both the structural details on very small scales and the 
associated very fast dynamics by using a coarse-graining along the chemical sequence of 
the chains [O]: n « 3 — 5 successive chemical units are mapped into one effective bond, 
connecting two effective monomers. These effective bonds may be defined on a lattice as 
well as in the continuum. In the lattice case, the effective monomers are formed from all 
8 sites of an elementary cube on the lattice, and no site may belong to more than one 
monomer (excluded volume interaction). Bond lengths in this "bond fluctuation model" 
[|TTH are taken in the range from 2 to \/l0 lattice spacings, and bond vectors b = (±3, 0, 0) 
[or permutations thereof] are taken to be in the ground state, while all other choices of 
the bond vector represent excited states, which an energy cost e(— 1; also ke = 1). This 
bond length potential may be thought of as representing the original atomistic intrachain 
potentials in a coarse-grained way. Since this particular choice of the potential has the effect 
that each bond that is in the ground state blocks the 4 sites between its adjoining effective 
monomers from further occupation, it has the effect that free volume in the dense system is 
"wasted" and a "geometric frustration" is created - at high densities not enough free volume 



is available that all bonds manage to get into their ground state as T — > [12]. This model 



can be simulated very efficiently, and for chain lengths N as small as = 10 properties 
typical of polymers (e g. gyration radius (R^) oc y/N) are reproduced. A density of 



<p = 0.533 of occupied sites corresponds to a dense melt |T^,[T3[. 

Dynamics is introduced by randomly selecting a monomer of a chain and by attempting 
to move it by one lattice unit in a randomly chosen direction. Of course, the move is only 
carried out if the new position satisfies excluded volume and bond length constraints, and the 
possible energy change is considered with the usual Metropolis criterion [LT],^2]|. It is thought 



that this "random hopping" of effective monomers in a coarse-grained sense corresponds to 
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the random hops of chemical bonds over barriers in the torsional potential, because by such 
moves pieces of a chain can rotate relative to each other, and thus are responsible for the 
relaxation of chain configurations |IT|,P1. If one carries out such a mapping of a realistic 



(i.e., atomistic) chain model on the bond fluctuation model literally, one estimates that one 
attempted Monte Carlo step (MCS) per monomer corresponds to 1CT 13 sec, and the lattice 
unit to a distance of 2A |14]]. Typically, runs are carried out for systems of L x L x L with 
periodic boundary conditions and L = 30, averaging over 16-160 independent "replicas" of 
the system over a time of up to 10 7 MCS. Thus even for the time window of this coarse- 
grained model we are restricted to t ~ 1/zsec (while molecular dynamics simulations of 
atomistic models can deal with t<lnsec only). However, one important aspect is that one can 
equilibrate the system configurations also with artificial moves, which have no counterpart 
in the dynamics of real chains, such as the "slithering snake" algorithm [OJTJJ, and thus 



gain several orders of magnitude in time due to faster relaxation. In this way, configurations 
of the well equilibrated melts are prepared at rather low temperatures, which serve as initial 
states for runs with the "random hopping" algorithm studying the dynamics of these melts. 

Alternatively, we treat a continuum model where an effective bond is represented by a 
(finitely extensible) spring, and the effective monomers are the beads of this bead-spring 
model. There is a Lennard- Jones interaction between any pair of beads cut off at a distance 



r c = 2 x 2 1 /V |5|, 



U u {r)=Ae (a/r) 12 -(a/r) 



+ c, 



(1) 



the constant c being chosen such that Uhj{ r c) — 0. Choosing here units of length and 
temperature such that a — 1, e — 1, the springs are described by 



U(r) 



2 



1 - (r/Ro) 



R = 1.5, k = 30. 



(2) 



While ULj(r) has its preferred distance at a minimum for r min = 2 1 / 6 , the minimum of the 
potential between neighboring monomers along the chain (given by the sum of Eqs. ([[]), (|2|)) 
occurs for tq « 0.96: although the bead-spring model is fully flexible, this incompatibility of 
fo/ r min with any standard crystal structure prevents crystallization of our model, and thus 
creates the "frustration" in the fluid necessary for the formation of a glass upon cooling. 

This off-lattice model is simulated by NpT-molecular-dynamics, choosing also a chain 
length of N = 10, and altogether M = 1200 monomers in the system Note that 

the constant pressure ensemble [I6| is used for equilibration only - for a study of dynamic 
properties a clock is set to zero and runs are started at T =const (using a suitable Nose- 
Hoover thermostat |17| ) which produce a dynamical behavior practically indistinguishable 
from the microcanonical ensemble 
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Fig. 1 shows that one obtains the structure factor S(q) of the amorphous polymer melt 
in close similarity with the experiment ||. If one studies the volume v(T) per effective 
monomer in a slow cooling run, starting at p — 1 from a well-equilibrated configuration at 
T = 0.6 and lowering T every 500 000 MD time steps by 0.02, one finds a rather well-defined 
kink at T, 



g ~ 0.41, while a fit of the self-diffusion constant to the Vogel-Fulcher law yields 



- VF 



o.33 m 



Thus the model does exhibit a glass transition, as expected. 
The bond fluctuation model also yields S(q) ~ for small q due to the very small 
compressibility of the polymer melt, and then a broad peak similar to the "amorphous 
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halo" of Fig. 1 occurs, at q ~ 3 (in units of the lattice spacing for this model; physically this 
corresponds to about 1.5 A ) JL3|. The peak position and shape also change very little with 
temperature. However, the second shallow and temperature-independent peak of Fig. 1 at 
still larger q, representing intra-chain correlations, is less well reproduced jOJ. This is not 
surprising, since q = 2ix corresponds in real space to one lattice unit. A further disadvantage 
of the lattice model is that one can study the glass transition only at constant volume, while 
in the off-lattice model both constant volume and constant pressure studies were performed 



III. TESTING THE MODE COUPLING THEORY (MCT) OF THE GLASS 

TRANSITION 

Fig. 2 shows the intermediate incoherent dynamic structure factor F q (t), defined as 

^ M 

= m ? (eXpW ' [m ~ m]}) ' (3) 

i=l 

where the sum is over all M monomers in the system (which are at position r*j(t) at 
time if:). Pronounced two-step relaxation is seen, the second step ( "a-relaxation" ) satisfies 
the time-temperature-superposition principle, while the first step ("beta relaxation") does 



not 1 15]. The relaxation time r q is both compatible with the Vogel-Fulcher law with Typ 



0.33, and with a power law (Fig. 3) but with a q-dependent exponent 7 g [fL8~l . While the 



consistency with a power law is evidence in favor of MCT @, the g-dependence of the 
exponent 7 g is not. In addition, one sees that the data deviate from the power law both for 
T/T c — 1 > 0.4 (as expected, since one has left the asymptotic region) and for T/T c — 1 < 0.04 
(which can be attributed to "hopping processes" by which effective monomers can escape 
from the cage formed by their neighbors, and requires use of the extended version of the 
theory || that describes the rounding of the ergodic to nonergodic transition). Thus, the 
observability of the idealized MCT M is restricted to about one decade in T/T c — 1, showing 



the limitations of this theory when applied to polymers ||15] , |I8| , |19[ . 

In the a- regime F q (t) can be fitted to the KWW law with an exponent (3 that is also 
weakly g-dependent {(3 q = 0.70 ± 0.08 for q = 6.9 [0); for a detailed analysis of this q- 
dependence see Ref. |19|). The predictions of MCT for the (3 regime can also be tested 
in detail and compare rather favorably with MCT p^| , although the same caveat over the 
restricted "temperature window" where MCT is applicable must again be made. 

This analysis can be replaced for several pressures and thus one can trace out a "critical 
line" T c (p) in the (T, p) plane separating the liquid from the ideal glass (Fig. 4). Carrying 
out simulations at constant pressure and at constant density that lead to the same point 
{T c (p),p) on the line, one finds that indeed r q oc (T — T c )~ 79 holds with the same T c and the 
same 7 g for both paths in the (T, p) plane. This underscores that the MCT critical line is 
indeed physically significant. 



Also the bond fluctuation model has been compared to MCT, both in its idealized fl20 



and extended Ell version. While a fit of various relaxation times and of the selfdiffusion 



constant to the VF law is nicely compatible with the data and yields Tvf ~ 0.125±0.005 [22 



the MCT fits both yield T c « 0.15, and power law behavior occurs over a similar 
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temperature interval as for the continuum model. One characteristic difference between 
both models concerns the /3-relaxation regime, however: while in the off-lattice model the 
first decay of the structure factor F g (t) {Fig. 2} is due to small-amplitude motions (over 
distances of the order of 10% of the distances between effective monomers and their nearest 
neighbors), no such small scale motion is possible in the bond fluctuation model: either a 
monomer hops a lattice unit), or it cannot hop at all. This distinction also shows up when 
we compare the time dependence of the scaled mean square displacements of monomers 
for both models (Fig. 5). Both the lattice model and the continuum model are thought 
to correspond essentially to the same physical system, a dense melt of short polymers, on 
a coarse-grained level. This implies that on mesoscopic scales (length scales of the order 
of the distance between effective monomers or larger) these different coarse-grained models 
should yield very similar results. Comparing suitably scaled data, so that units of length 
and time are absorbed in the scales, this is indeed the case (Fig. 3). Note that the regime 
where MCT is applicable is actually the regime of rather small displacements, up to the 
order of distances between monomers, from where on then a Rouse-like behavior controls the 
dynamics. 

IV. COMMENTS ON OTHER THEORETICAL CONCEPTS AND SOME 

CONCLUDING REMARKS 

A particular advantage of the lattice model is that the slithering snake algorithm allows 
to obtain equilibrium properties of the model at fairly low temperatures, and also the con- 
figurational entropy S(T) can be computed [^] in order to test the Gibbs-Di Marzio theory 
||. Indeed one finds that S(T) decreases to about 1/3 of its high temperature value S(oo) 
as T approaches T g , but there the curve S(T) vs T bends over and gives evidence that S(T) 
stays nonzero at lower temperatures. Conversely, if one extracts the quantities that are used 
in the theory || directly from the simulation and inserts them into the (approximate!) the- 
oretical formula one would find an "entropy catastrophe" {S(T) becoming negative} 
at T < 0.18. This unphysical result, however, is easily traced back to a severe underesti- 
mation of S(oo) by the approximations of ||. The "entropy catastrophe" |J thus clearly is 
an artefact of an inaccurate approximation, one should not attribute physical significance 
to it. In such manner simulations can go beyond experiment for testing theories, because 
the input parameters of the theories can be unambiguously extracted from the simulation 
as well, and a comparison between theory and simulation is possible without adjustable 
parameters whatsoever. 

It is interesting to note that nevertheless the Adam-Gibbs equation []7), D/D oc 
exp{— i^act /TS(T)}, where D is the self-diffusion constant at very high temperatures and 
E'act some activation energy, provides a very good description of the diffusion constants found 
in the simulation when one uses also the S(T) found in the simulation (Fig. 6). At the same 
time, the absence of any finite size effect in Fig. 6 is at odds with the idea to attribute the 
slowing down implied by the decrease of D(T) to an increasing glass correlation length £(T) 
0, via r(T) oc [£(T)] 2 with z a dynamic exponent fllOl : if such an hypothesis would hold, one 
would expect that decreasing L should decrease £(T) [since £(T) < L] and hence r(T), and 
in turn D(T) should increase with increasing L. Fig 6 proves the absence of such finite size 
effects. This finding is surprising, since a growing static length can be extracted both from 
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the pair correlation function in the melt |24| and from surface effects near hard walls |]25| 



It hence appears that there is a growing length identified in p^ , |25[| but it is not responsible 
for the slowing down at T g . Of course, there is no contradiction with experimental results 
which find that T g in thin films (or in pores, respectively) changes when the linear dimension 
of the film (or pore respectively) is varied: depending on the boundary conditions at the 
surface, the local mobility of monomers in the surface region changes, and this effect is the 
more pronounced on the freezing the smaller the linear dimension. No such surface effect is 
present with the periodic boundary conditions in Fig. 6, of course. In experiments on the 
glass transition in confined systems, finite size effects and surface effects can never be clearly 
separated, unlike simulations where one can show that there is no finite size effect (Fig. 6) 
at least in the temperature region studied here, although there do exist surface effects . 

Returning to MCT, we emphasize that idealized MCT does provide a good description 
of a large number of simulation data, but only over a rather restricted range of temperatures 
(about one decade in T/T c — 1) and corresponding times (or viscosities, typically the range 
10° < rj(T) < 10 2 Poise if r](T c ) = 10 3 Poise: the huge range from 10 3 < r] < 10 13 Poise is 
outside the scope of the theory). Correspondingly, only a small intermediate regime of small 
monomer displacements in Fig. 5 is described - neither the initial increase that depends 
on microscopic properties of the model, nor the regime of hopping processes that lead to 
a Rouse-like relaxation of coil configurations (before ultimately ordinary diffusion sets in) 
are part of the theory. Since there occur smooth crossovers between the various regimes, a 
reliable assessment of the validity of MCT for polymer melts near T g is difficult. A more 
complete theory (that unifies e.g. MCT and the Rouse model) would be very desirable. 
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FIGURES 



FIG. 1. Structure factor S(q) plotted vs wavenumber q, for a system of 120 off-lattice 
bead-spring chains with chain length N = 10, simulated in the NpT-ensemble at scaled pres- 
sure p = 1, choosing Lennard- Jones units e = l,a = 1 (and /cb=1). Note that the zero of the 
ordinate for each curve is shifted upward by 0.001 relative to the previous one. From Bennemann 
et al. p|. 



FIG. 2. Dynamic structure factor F q (t) of the off-lattice model at scaled pressure p = 1 plotted 
vs a rescaled time t/r q (where T q is defined from F g (t = r g ) = 0.3) for q = 6.9, the peak position 
of the static structure factor at low T (cf. Fig. 1). Only temperatures in the range 0.46 <T< 0.7 
are included, as indicated. From Bennemann et al. fl5|| . 



FIG. 3. Log-log plot of the relaxation time r q versus T — T c with T c = 0.45, for a pressure p = 1 
and three values of q. The straight lines indicate the power laws T q oc (T — T c )~ lq with exponents 
7 g = 2.3, 2.1 and 2.0 (from above to below). From Bennemann et al. flj~8"| . 



FIG. 4. Critical line T c {p) in the (T,p) plane, for the off-lattice bead spring model. From 



Bennemann et al. 18 



FIG. 5. Comparison of the mean-square displacement gi(t) of inner monomers 
{9i{t) = ([n(t) -fj(0)] 2 J)}plottedvs scaled time {Dt/(R 2 g }, D being the selfidiffusion constant and 
{Rg) the mean square gyration radius of the chains} for the off-lattice model (curve marked MD) 
with the corresponding Monte Carlo results for the bond fluctuation model at three temperatures. 



FIG. 6. Selfdiffusion constant of the bond fluctuation model plotted vs inverse temperature, 
for L x L x L lattices with periodic boundary conditions and several lattice sizesL. The solid and 
the dotted line are fits to the Vogel-Fulcher and Adam-Gibbs equation, respectively. Data taken 
from Binder et al. 1241. 
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